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1 Introduction 

In many simulations of astrophysical objects and phenomena, the equation of state 
(EoS) of dense matter is an essential ingredient. It determines, e.g., the dynamical 
evolution of core-collapse supernovae BE! and neutron star mergers [3j, and the 
structure of compact stars [4]. The application of an EoS is reasonable if the timescales 
of reactions are much smaller than those of the system evolution and thermodynamic 
equilibrium can be assumed to hold. In general, a global EoS is required that covers 
a wide range in temperature, density and isospin asymmetry. These conditions affect 
the chemical composition of matter and the nucleosynthesis. 

A critical examination of existing global EoS models [5] suggests that the devel¬ 
opment of an improved EoS is worthwhile. The set of constituent particles should be 
enlarged considerably including not only nucleons, charged leptons and photons but 
also a “complete” table of nuclei, mesons, hyperons or even quarks as degrees of free¬ 
dom at high densities and temperatures. The model parameters have to be contrained 
better taking, e.g., properties of nuclei, results of heavy-ion collisions or compact star 
observations into account. Correlations should be considered more seriously, e.g., 
at low-densities where the virial equation of state (VEoS), which is determined by 
nucleon-nucleon correlations, is a model-independent benchmark 0IZI- For compos¬ 
ite particles such as nuclei the dissolution in the medium (Mott effect) has to be 
described properly m. Electromagnetic correlations are essential in order to model 
the solidification/melting at low temperatures. Phase transitions and the appear¬ 
ance of ’non-congruent’ features have to be treated correctly ra with non-negligible 
differences between nuclear matter and stellar matter. Obviously, it is a tremendous 
challenge to cover the full range of thermodynamic variables in a single unified model. 

Information on correlations are encoded in spectral functions, which have a com¬ 
plicated structure in general. Often, a quasiparticles (QP) approach is employed as an 
approximation. The QP properties change inside the medium and the size of residual 
correlations is reduced. The QP concept is very successful in nuclear physics, e.g., in 
phenomenological mean-field models (Skyrrne, Gogny, relativistic) or the treatment of 
pairing correlations using a Bogoliubov transformation m ■ In the ultimate limit, an 
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exact diagonalisation of the Hamiltonian of the interacting many-body system leads 
to a system of independent QP that can be many-body states. At low densities, clus¬ 
ters appear as new degrees of freedom as described in the VEoS. In order to consider 
these features, a generalized relativistic density functional (gRDF) was developed. It 
takes the correct limits and explicit cluster degrees of freedom into account. 

2 Generalized relativistic density functional 

The gRDF model [12, [[3] E] H5] is based on a grand canonical approach. It is 
an extension of a conventional relativistic mean-held model with density dependent 
couplings [16]. All thermodynamic quantities are derived from a grand canonical 
potential density co(T, {/q}), which depends on the temperature T and the set of 
chemical potentials /q of all particles. The present set of particle species comprises 
baryons (nucleons and hyperons), nuclei, charged leptons and photons. Besides light 
nuclei ( 2 H, 3 H, 3 He, 4 He) a full table of heavy nuclei ( A Z with A > 4, N, Z < 184) is 
included, too. Experimental binding energies are used or, if not available, predictions 
from the DZ model na. Internal excitations of heavy nuclei are considered with 
temperature dependent degeneracy factors obtained with appropriate level densities. 
Effective continuum resonances represent nucleon-nucleon scattering correlations and 
ensure the correct low density limit in accordance with the VEoS ra- 

All massive particles are treated as QP with scalar (S/) and vector (R) poten¬ 
tials. The effective interaction is modeled by an exchange of mesons (cr, uj, p) with 
density-dependent couplings to the nucleons, both free and bound in nuclei, using the 
well constrained DD2 parametrization P2J. It gives very reasonable nuclear matter 
parameters at a saturation density of n sat = 0.149 frn -3 , such as a binding energy per 
nucleon E/A = 16.02 MeV, a compressibility K = 242.7 MeV, a symmetry energy 
J = 31.67 MeV and a slope parameter L = 55.04 MeV. The neutron matter EoS lies 
within the error bounds of recent chiral effective held theoretical calculations [18] [19] . 
Both potentials S) and V] receive contributions from the meson fields. For composite 
particles, the scalar potential contains an additional mass shift Am, that depends 
on all particle densities and temperature. It mainly takes the blocking of states 
by the Pauli exclusion principle into account and serves to describe the dissolution 
of clusters by reducing the particle binding energy. This microscopically motivated 
approach replaces the traditional, purely geometric concept of the excluded-volume 
mechanism m- The vector potential V] includes a “rearrangment” contribution due 
to the density dependence of the meson-nucleon couplings, which is required for the 
thermodynamic consistency of the model, and an electromagnetic correction to ac¬ 
count for electron sceening effects in stellar matter. 
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3 Symmetry energy and neutron skins of nuclei 

The isospin dependence of the effective interaction in the gRDF model determines 
the density dependence of the symmetry energy. It is crucial for a proper description 
of the structure of neutron stars, see, e.g., the topical issue on the symmetry energy 
[21] . A strong correlation of the neutron skin thickness A r np of heavy nuclei with 
the slope of the neutron matter equation of state [22, [23] or the slope parameter 
L of the symmetry energy is observed when the predictions of a large number of 
mean-held calculations, both relativistic and non-relativistic, are compared, see, e.g., 
[21J. In recent years, many attempts were made to determine the symmetry energy at 
saturation J and the parameter L from experiments, e.g., by measuring the neutron 
skin thickness of 208 Pb and using the A r np vs. L correlation. Since the calculations of 
neutron skin thicknesses are based on mean-held models, the question arises whether 
few-nucleon correlations can effect the results. 

The gRDF approach can be employed to describe the formation of nuclei inside 
matter at hnite temperatures by using an extended Thomas-Fermi approximation in 
spherical Wigner-Seitz cells [13]. In a calculation with nucleons and light clusters 
as degrees of freedom, it is observed that the probability of finding light clusters is 
enhanced at the surface of the heavy nucleus as compared to the surrounding low- 
density gas. The gRDF model can be extended to the description of heavy nuclei 
in vacuum at zero temperature to study cluster correlations. In this case, only the 
a-particle remains as the relevant light cluster. Its density distribution is obtained 
from the a-particle ground state wave function that is calculated self-consistently 
in the WKB approximation. For the chain of Sn nuclei, a distinct reduction of the 
neutron skin thickness is observed when a-particle correlations are considered ra- 
However, the effect vanishes for very neutron-rich nuclei or for nuclei with roughly 
the same neutron and proton numbers without a neutron skin. A variation of the 
isovector dependent part of the effective interaction allows to study the A r np vs. L 
correlation, e.g., for a 208 Pb nucleus. A systematic shift is observed that might affect 
the determination of the slope parameter L from measurements of the neutron skin 
thickness, at least as a systematic error. It is envisaged to investigate experimentally 
the predicted formation of a-particles at the surface of Sn nuclei in quasi-elastic (p,pa) 
reactions at RCNP, Osaka [26]. 


4 Outlook 

The present version of the gRDF model includes only hadronic and leptonic degrees 
of freedom where nuclei are described as clusters composed of nucleons. At high den¬ 
sities or temperatures a phase transition to quark matter is expected. Hence, quark 
degrees of freedom should be incorporated into the approach. On the other hand, at 
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low densities and temperatures, quarks should be confined in nucleons. In a prelim¬ 
inary extension of the gRDF model with quarks, a phenomenological description of 
confinement will be implemented. The idea is to apply an “inverse” excluded-volume 
approach that permits the quarks to propagate freely only above a certain (scalar) 
density of the system. For this purpose, the classical excluded-volume mechanism is 
generalized by allowing more general dependencies of the “available volume fraction”. 
The correct quantum statistics and a relativistic description are considered, too. The 
relevant theoretical formulation to guarantee the thermodynamic consistency of the 
approach has been developed and exploratory calculations have been preformed. 

Another extension of the gRDF model concerns the introduction of more general 
meson-nucleon couplings in the Lagrangian density. In conventional RMF approaches 
with density-dependent couplings, the nucleon self-energies only depend on densities. 
As known from Dirac-Brueckner calculations of nuclear matter, they should also de¬ 
pend on the nucleon momentum or energy. This dependence can be mapped to 
modified effective density dependent meson-nucleon couplings [27], but the full de¬ 
pendence should be kept in order to comply with the optical potential constraint at 
high nucleon energies. This can be achieved in a RMF model with density-dependent 
and non-linear derivative meson-nucleon couplings of general functional form [28] . 
Preliminary studies indicate a softening of the EoS at high densities, however, for a 
reliable fit of the model parameters, the approach has to be applied to the description 
of finite nuclei. Work in this direction is in progress. 
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